From time inversion to nonlinear QED* 



Wei Min Jin 

Department of Physics and Astronomy, State University of New York at Buffalo, 
Buffalo, NY 14260-1500, U.S.A.^ 
Electronic mail: wjinl@excite.com 



In Minkowski flat space-time, it is perceived that time inversion is unitary rather than an- 
tiunitary, with energy being a time vector changing sign under time inversion. The Dirac 
equation, in the case of electromagnetic interaction, is not invariant under unitary time 
inversion, giving rise to a "Klein paradox" . To render unitary time inversion invariance, a 
nonlinear wave equation is constructed, in which the "Klein paradox" disappears. In the 
case of Coulomb interaction, the revised nonlinear equation can be linearized to give en- 
ergy solutions for Hydrogen-like ions without singularity when nuclear number Z > 137, 
showing a reversed energy order pending for experimental tests such as Zeeman effects. In 
non-relativistic limit, this nonlinear equation reduces to nonlinear schrodinger equation 
with soliton-like solutions. Moreover, particle conjugation and electron-proton scattering 
with a nonsingular current-potential interaction are discussed. Finally the explicit form 
of gauge function is found, the uniqueness of Lorentz gauge is proven and the Lagrangian 
density of quantum electrodynamics (QED) is revised as well. The implementation of 
unitary time inversion leads to the ultimate derivation of nonlinear QED. 
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1. INTRODUCTION 



Decades ago, Wigner [1] first pointed out the significance of time inversion and later 
made detailed discussions on antiunitary time inversion (Wigner [2]). His theory about 
time inversion is based on a classical motion picture (Wigner [3]): "Time inversion . . . 
replaces every velocity by the opposite velocity, so that the position of particles at +t 
becomes the same as it was, without time inversion at — t. ..." 

Also well known is Einstein's relativity theory (Einstein [4]) that completely changes 
our perception of space and time, as best manifested by Minkowski's "world-postulate" 
(see Minkowski's paper in [4]). Minkowski unifies space and time by introducing time as 
an independent coordinate in addition to three space coordinates. Usually "Minkowski 
space-time" refers to four-dimensional flat space-time in special relativity, in which rela- 
tivistic quantum mechanics and quantum electrodynamics (QED) are established. 

On the development of quantum mechanics and relativity, there always exist heated 
debates on the physical meaning of the theories. Dirac however was indifferent in such 
debates. Instead he had been looking for mathematical possibilities to reconcile quantum 
mechanics and relativity. As a consequence, he founded relativistic quantum mechanics 
with a relativistic wave equation called "Dirac equation" (Dirac [5]). However, the Dirac 
equation was also not perfect. In the following year, Klein [6] found a paradox in the 
Dirac equation. Now seven decades has gone by, the "Klein paradox" has still not been 
resolved mathematically, though it can be explained away by physical reasoning. 

This paper is organized as follows. Given that a correct concept in Newtonian classi- 
cal mechanics may not be necessarily correct in Einstcinian relativistic mechanics, first I 
am going to show that the above mentioned classical motion picture is just one of those 
concepts, correct classically but not relativistically. Then I will proceed to clarify what 
is supposed to mean to make a time inversion in special relativity. Based on the general 
principles of quantum mechanics and special relativity, it can be shown that time inver- 
sion turns out to be a unitary transformation with energy being a time vector changing 
sign under time inversion in Minkowski space-time. 

With unitary time inversion understood, the next step is to study the transformation 
of the Dirac equation. For clarity, let us name the Dirac equation in the case of electro- 
magnetic interaction as Dirac EM-equation to distinguish it from the Dirac equation in 
other cases. What happens is: the Dirac EM-equation is not invariant under unitary time 
inversion. As a direct consequence of this non-invariance, the Dirac EM-equation has a 
non-symmetric positive-negative energy spectrum with a mathematical singularity that 
causes a "Klein paradox" . Here in this paper, I do not intend to argue too much about 
whether or not the "Klein paradox" has physical meaning. Rather I hold such a point of 
view that we would be better off from the outset without mathematical singularity. 

To implement unitary time inversion, I come up with a revised equation that looks 
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similar to the Dirac EM-equation on the one hand, while differs from the Dirac EM- 
equation in several aspects on the other. The main difference is: the revised equation 
preserves the invariance under unitary time inversion in addition to the invariance under 
the other Lorentz transformations. Consequently, it gives a positive-negative symmetric 
energy spectrum without singular crossing point. In the case of Coulomb interaction, the 
energy bound states for Hydrogen-like ions are solved without singularity when nuclear 
number Z > 137, and the order of some energy levels found opposite to the conventional 
one, which entails further high precision tests from experiments like Zecman effects. 

Another difference is; the Dirac EM-equation, involving the external electromagnetic 
potential with minimal coupling, is linear on fermion field; while the revised wave equa- 
tion, involving the electro-dynamical interaction potential, is nonlinear on fermion field. 
In non-relativistic limit, the Dirac EM-equation reduces to the linear Schrodinger equa- 
tion, while the revised nonlinear wave equation reduces to the nonlinear Schrodinger 
equation with soliton-like solutions well known in nonlinear physics. 

Typically, when the finite size of fermion is considered, the interaction potential can 
be generalized by a convolution between four-current and four-potential with an integral 
over the finite four-dimensional size of the fermion. Using such a nonsingular convolution 
can avoid the singularity and infinity troubles in calculating the self-energy of electron 
and the divergent integrals of Feynman Diagrams. 

The gauge invariance is to be discussed at the end. To preserve Lorentz invariance 
including the invariance under unitary time inversion, it can be shown that gauge func- 
tion is of particular form, not arbitrary, and the only gauge condition is Lorentz gauge. 
It is also straightforward to write a nonlinear Lagrangian that derives the Maxwell elec- 
tromagnetic field equation and the revised nonlinear fermion field equation. With the 
basic principles of quantum field theory, we may then establish a nonlinear QED that 
would show many interesting applications in experiment. What I am actually up to, is to 
make necessary improvements on QED in its own framework, namely in the language of 
space-time, without further drastic changes as in superstring theory. So much has been 
said, let me now turn into more detailed discussions. 



Usually under a coordinate transformation x' = ax, the transformation of wave function 
in quantum mechanics is defined by 



2. UNITARY TIME INVERSION 



A(a)*(a;) = *(aa;). 



(2.1) 



But antiunitary time inversion is defined by (see Bjorken and Drell [7]) 



T**(t,x) = **(-t,x). 



(2.2) 
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where the extra complex conjugation * on wave function is used to comply with the 
assumption that under time inversion (t.x) (— f.x), energy and momentum transform 
as {E, p) — > {E, — p). In this case, the phase of a plane wave (p = p - x — Et (natural 
units are used in this paper), is not invariant under antiunitary time inversion. This also 
leads to another result that the imaginary unit i has to be supposed to change to — i 
under antiunitary time inversion, while constant i has nothing to do with time. 

In quantum mechanics, the group velocity of a wave packet is expressed by dx/dt = 
dE{p)/dp. By making a positive non-relativistic energy assumption E{p) = p^/2m, 
one can then get a classical correspondence dx/dt = p/m = v, from which the classical 
motion picture comes. The point is: in the derivation of classical mechanics velocity 
from quantum mechanics group velocity, the positive energy assumption has been taken, 
which may not be correct under time inversion in special relativity. It would be logically 
odd to discTiss symmetry under time inversion based on a truncated classical motion 
picture that has already broken the symmetry. 

Now that quantum mechanics is more general than classical mechanics, it is better 
not to impose any properties for time inversion classically until we obtain certain results 
from quantum mechanics. As we understand, time inversion is nothing but a kind of 
coordinate transformation in space-time, and it is natural to adapt the general definition 
(2.1) for time inversion rather than to follow the truncated classical motion picture. 
Therefore putting the general principles of quantum mechanics in the first place, we 
define time inversion as; 

T'i'(t,x) = *(-f,x). (2.3) 

For a plane wave \E'(t, x) = C(p, E) exp[z(p • x — Et)], it is clear that the expectation 
value of space momentum operator p = — i9x remains the same and that of energy 
operator E = idt changes sign under time inversion. Any wave function can be Fourier- 
transformed into a combination of plane waves. So it is easy to check these results hold 
for any kind of matter waves. 

In 4-d space-time, the space-time interval squared AX^AX^ ~ g^'^AX^AX^ is 
Lorentz invariant, where 5°° = —g" = l{i = 1,2,3) and g'^'^ = ^ v). Mathe- 
matically if contravariant four-coordinate transforms as X'^ = a^^X", then covariant 
four-coordinate transforms in the way X'^ = aJ^X^. Since AX'^AX'^ = AX^AJf , we 
get = ga'^g where a~ is the transpose of a, leading to det(a~^) = det(a) = ±1 
for proper and improper Lorentz transformations respectively. On the other hand, the 
phase of a plane wave (j) = —P^X^ characterizing the physical correlation of space-time 
world-points, is also invariant under homogeneous Lorentz transformations, therefore 
four-momentum must transform in a covariant way = a'/^Pv In the case of time 
inversion (t, x) (— t, x), it is naturally found an energy-momentum transformation 

{E,p)^{-E,p). (2.4) 
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Moreover, if a plane wave with a constant four-momentum = (E, p) travels from 
space-time point 1 to 2, the phase difference between them Acp = 02 — <Pi describing the 
causality of this process, is an invariant under all inhomogeneous Lorentz transformations 
X'^ = a^X'^ + b^. The phase differences between points 1 and 2 before and after 
space-time inversions can be written down as: = p- (x2 — Xi) — E{t2 — ti) and 
A^' = p' • (xj — x'j^) — E'{t'2 — t[). Since = A(j), we obtain p' = — p under space 
inversion xi^ — x'j^ = — (x2 — Xi), and E' = —E imder time inversion ^2 ~ *i = ~(^2 — ii)- 
It is such a Lorentz invariant physical phase space that determines the unique energy- 
momentum transformation (-E,p) {—E,p) under time inversion. 

Actually four- vector momentum = {E, p) consists of time component E and space 
component p. While momentum p is usually known as a space vector changing sign under 
space inversion, energy E may be similarly regarded as a time vector changing sign under 
time inversion. The idea of uncertainty principles in quantum mechanics: ApAx ~ 1 
and AEAt 1, is that the strict measurement of momentum and position as well as that 
of energy and time can not be taken simultaneously since they do not commute, but they 
do have close relationship with each other. The determination of momentum is tightly 
related to that of position not to that of time. Similarly the measurement of energy is 
directly connected to that of time not to that of position. This one-to-one correspondence 
reveals the following fact that 4-d space-time and energy-momentum, combined into a 
Lorentz invariant phase space, are the reciprocal systems relatively to each other 

= id^,, (2.5) 

with a commutator relation [X'',P^] = —i. 

The Fourier transformation used in quantum mechanics also shows a good example of 
this one-to-one correspondence. According to this correspondence, we prefer to call 4-d 
energy-momentum system, the "reciprocal world" of space-time. Generally speaking, this 
is only a question of using different representations or interchangeable languages. There 
is nothing special, in the sense of relativistic covariance, that producing an inversion 
in space-time world is equivalent to making a corresponding inversion in its energy- 
momentum reciprocal world by the hint of (2.5). This idea now enables us to understand 
why momentum switches sign under space inversion while energy reverses sign under 
time inversion. 

The physical picture becomes clear if in a "local framework", we do not impose 
absolute future and past concepts by setting a standard clock A but rather reverse time 
axis by setting another clock B running counterclockwise. If a plane wave is propagating 
toward the future with a positive (or negative) frequency by the standard clock A, then 
it can also be equivalently looked upon as propagating toward the past with a negative 
(or positive) frequency by the other clock B. On the principle of special relativity, it is 
arbitrary in setting clocks in locally flat space-time. Therefore this plane wave may have 
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either positive or negative energy depending on which clock we use. 

By the correspondence between space-time and energy-momentum worlds, if we in- 
deed want to introduce a time inversion concept in Lorentz group, then it is natural to in- 
troduce a corresponding "energy inversion" concept to map these two worlds completely. 
Once again we would like to emphasize, based on what we try to clarify here from the in- 
trinsic structures of space-time world and energy-momentum reciprocal world, and from 
the general principles of special relativity and quantum mechanics, that in Minkowski 
flat space-time, energy is not a scalar any more, it is a "time vector" changing sign under 
time inversion. 

Similar to space inversion, by definition (2.3) it is easy to prove: (a) T is linear; (b) 
= 1; (c) T* =T; (d) T-" =T. So time inversion operator is unitary 

= T = (2.6) 

To find its eigenvalues let = A^*. Prom relation (b) we get A = ±1 representing even 

and odd "time parities" respectively. 

For a free relativistic particle with an energy- momentum relation = 
the common eigenfunctions of momentum operator p = —id^ and energy operator E = 
idt are expressed by $±(t,x) = C{p,E)(j)p{x)exp{±iEt). Obviously they are not the 
eigenfunctions of time inversion operator T since T does not commute with E: [T, E] ^ 0. 
Now make linear combinations: ^'-i- = (<!>+ + $-)/2 and = ($_|_ — $_)/2i, which are 
the common eigenstates of p and T but not those of E any more. Generally any state 
wave function can be divided by ^1/ = + ^f*-, where the even and odd time parity 
eigenstates are 'i!± = (1 ± T)\l//2 respectively. And any operator can be expressed by 
W = W+ W_, where W± = (W ± TWT)/2 and TW±T = ±W±. Here W+ represents 
the even time parity operator such as momentum p = —id^ and W_ represents the odd 
time parity operator such as energy E = idf. 

3. DIRAC EQUATION REVISITED 

At the end of 19th century, Zeeman [8] detected the splitting of spectral lines under 
the influence of external magnetic fields, revealing the existence of electron spin. To 
explain the Zeeman effect, the imagination of electron structure has been made: electron 
has a magnetic moment and thus has an intrinsic spin angular momentum instead of a 
classical point particle (Lorentz [9]). It has been clarified that electron spin relates to 
the extra degrees of freedom and can be well described by a two component complex 
variable called "spinor". One is then faced with such a typical question: how to con- 
struct a minimal complete set of variables in the combination of space-time and intrinsic 
spin. In Dirac's mind (Dirac [10]), spinors, like tensors, are geometrical objects, yielding 
covariant transformations with respect to Minkowski space-time, and specifically under 
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Lorentz transformations, the Dirac equation is to obey the covariance principle of special 
relativity. Along this line of thought, we may say the whole domain of space-time and 
intrinsic spin is Lorentz invariant, and name it "common space" for clarity. 
A common variable is now defined by 

^ = YX^, (3.1) 

where matrices 7''(/U = 0,1,2,3), in the representation of Bjorken and Drell [7], are 
ascribed to the spin geometrical factors. The common variable constructed in this way is 

a function of space-time coordinates and also a metrical quantity. The common variable 
interval squared becomes AO^ — AX^AX^, leading to the anticommutation relations 
of Dirac matrices |7'',7''| = Sg'^''. Moreover, the derivative with respect to common 
variable can be deduced as 

d. = {-^)d^=l^d,, (3.2) 

and a momentum operator in common space is naturally introduced as 

= id^ = ji^P^. (3.3) 

By (3.1) and (3.3) it is straightforward to prove a commutator relation 

[n,P^] = -Ai. (3.4) 

If noticing that the common momentum squared of a free particle with mass m is constant 
P^ = — = m? , we find P^ has two eigenvalues ±m. If choosing +m eigenvalue, 
we arrive at the Dirac equation of free spin one-half particle 

^i^P^^{x) = m*(a;), (3.5) 

The Dirac equation, on the principle of special relativity, is required to be invariant 
under Lorentz transformation L(^ in common space, which is a direct product of spinor 
one Lg and coordinate one Lc: L^^ = L^Lc- Making on both sides of (3.5) and noting 
that Lg commutes with space-time vectors and Lc commutes with spinor vectors, we have 

L,7''L7iLcP^L-iL^*(a;) - mL^^{x). (3.6) 

Since momentum P^ is a covariant four-vector in space-time, 

LeP^L-i = a^P^, (3.7) 
matrices 7^ must be correspondingly contravariant in spinor space, 

L,7'*L;1 = 0(^7^ (3.8) 
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The above discussions also hold for two special cases. First, there is a unitary space 
inversion in common space S^; = SgSc, where = 7° is a unitary space inversion in 
spinor space: 

S,7°S;i = 7°, S.7SJ^ = -7, (3.9) 

and Sc is a unitary space inversion in space-time, defined as usual in quantum mechanics. 
Second, there is a unitary time inversion in common space = T^Tc, where Tg = 
7^7^7^ is a unitary time inversion in spinor space: 

T,70T;1 = -7°, T,7T;1 = 7, (3.10) 

and Tc is a unitary time inversion in space-time, defined by (2.3). 

When a spin one-half particle is put into an external Maxwell electromagnetic field, 
the Dirac EM-equation with minimal coupling can be written down in a general form 
(Dirac [5] [11]): 

7''(P^ - eA^)*(ar) = m^{x), (3.11) 

where = ($, A) is the space-time four-potential of electromagnetic field, proposed by 
Minkowski [12]. Under time inversion, time component $ does not change, while space 
vector A is supposed to change sign. The Minkowski four-potential constructed in this 
way is therefore "time-anticovariant" under time inversion: 

($,A)^($,-A). (3.12) 

If four-momentum is also "time-anticovariant" , one may make spinor space yield 
a sort of transformation to keep the equation invariant as done in conventional theory. 
However from what we have shown in the last section, four-momentum is not "time- 
anticovariant", but rather "time-covariant" (2.4), just as it is covariant under all the 
other Lorentz transformations. Here exists a nontrivial dilemma that the summation 
of "time-covariant" four-momentum and "time-anticovariant" Minkowski four-potential 
multiplied by charge and a negative sign is neither time-covariant nor time-anticovariant, 
keeping in no way the Dirac EM-equation invariant under time inversion. 

This dilemma gives rise to an unsymmetrical positive-negative energy spectrum when 
we solve this equation. Consequently when the electric potential $ such as Coulomb 
potential is strong enough, the positive and negative energy spectra will be shifted to 
overlap, leading to a mathematical singularity difficulty - "Klein paradox" (Klein [6]). 
Particularly, for a Hydrogen-like ion there is no real bound energy solution but appears 
a singularity when the nuclear number of the ion is larger than 137. 

At first, Dirac could hardly understand why his equation has negative energy solutions 
that seem "nonphysical" . To find an explanation, Dirac constructed a hole theory with 
a vacuum of fully filled negative energy states (Dirac [13]). The predicted positron 
was found a few years later (Anderson [14]). Despite Dirac's profound prediction of 
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antiparticles, the hole theory has difficulties such as infinite density, infinite negative 
energy and vacuum fluctuation in the vacuum state. Based on the hole theory, the "Klein 
paradox" is explained in terms of vacuum polarization accompanying particle-antiparticle 
pair production and annihilation under strong fields. Along this line of work, one has 
been trying to explain away the "Klein paradox" rather than to avoid it. Eventually, one 
is led to such a conclusion: " 'Klein paradox' is not a paradox" . 

Just as in classical mechanics, energy in the Dirac hole theory is still viewed as a 
scalar: namely negative energy is always "lower" than positive energy. If energy in 
special relativity is a frame-dependent "time vector" rather than just a scalar, as we have 
shown earlier, then it makes little sense to insert negative energy into vacuum or explain 
the "Klein paradox" as vacuum polarization. There should exist a more fundamental 
solution for these problems. In the next section, we would like to provide a down-to- 
earth approach, to remove the "Klein paradox" by revising the Dirac EM-equation. It 
sounds a little radical. Indeed this can not be done without the understanding that time 
inversion is a unitary transformation and energy is a "time vector" changing sign under 
time inversion in Minkowski space-time. 

4. NONLINEAR EQUATION RECONSTRUCTED 

As shown in §2, energy is a time vector changing sign under time inversion. Therefore 
the whole energy spectrum we obtain from whatever relativistic wave equation we might 
construct, must be positive-negative symmetric. This is, on the other hand, a necessary 
check if the constructed equation is truly invariant under unitary time inversion or not. 
By looking into the Dirac EM-equation, we find that the non-symmetric energy spectrum 
is caused by a matrix [3 in front of electric potential $. If /? is removed there, the equation 
becomes invariant under unitary time inversion. But it is no longer invariant under the 
other proper Lorentz transformations, if $ and A are still external potentials. After care- 
ful checking, we realize that interaction potentials have to be utilized in lieu of external 
potentials in a special way, to render the invariance under all Lorentz transformations. 
The following equation is what we get: 

(7'*P^ - + 67 • A^)*(a;) = m^ix), (4.1) 

where and A'' are the interaction potentials related to the external potentials ^ and 
A. We are going to derive their relations by considering Lorentz transformations, among 
which time inversion is defined in the unitary way we presented earlier. 

As in §3, one can see this modified equation is invariant under unitary space inversion 
Sui — 7"Sc and time inversion T^^ = 'y^j'^'j^Tc where Tc is defined by (2.3), by assuming 
that scalar potential does not change sign under space-time inversions while space 
vector potential A^ changes sign under space inversion but not under time inversion. 
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These assumptions can be directly verified after we derive the expUcit forms of and 
a' in terms of $ and A. 

Below we are going to prove that the modified equation (4.1) does give a symmetric 
positive-negative energy spectrum, to enhance our confidence that this equation is indeed 
invariant under unitary time inversion. Prom the equation (4.1), the energy operator is 
expressed in the way 

E = pi<T-(p-eA^)+P3(m + e$^), (4.2) 

where 

P^=\ T n)'/'2={,r ~n )>P3=(n (4-3) 



^ I j"^^ \ ii y "^-^ \0 -I ^ 

form a set of 4 X 4 matrices analogous to 2 x 2 Pauli matrices ai{i = 1, 2, 3). Here / is a 
2x2 unit matrix. To diagonalize (4.2), take a look at the square of E 

= (p - eA^)2 - e<T • (V X A' + A' X V) + (to + e^^)^ - epao" • [V, $^]. (4.4) 



Making a unitary transformation (U ^ =U^) 

V=±( ' 



(4.5) 



on both sides of (4.4) we can diagonalize it by using 

= -P3, (4.6) 

to the following form 

r=uEV^ = f^' 

= (p - eA^)2 - eor • (V X A^ + A' X V) + (to + e^^)^ + epsa ■ [V, $^]. (4.7) 
If supposing the following solutions 

Kxu = Elxu, e'x; = Efxi, (4.8) 

we can rewrite (4.7) as 

e" = ^(IXu >El< Xu\ + \Xi > Ef < XiD- (4.9) 

X 

Transforming back to 

i^u = V-\u, i^i = V-\i, (4.10) 

we get 

e' = ^(iV-u >El< + IVi > Ef < (4.11) 



10 



It is straightforward to check that the energy operator turns out to be 

E = J^d^u > XEu < i'u\ + > \Ei < Vil), (4.12) 

with A = ±1, which indeed gives a symmetric positive- negative energy spectrum. 
To keep (4.1) invariant under the proper Lorentz transformations 

-J-' = Lp*,*' = *L-\ (4.13) 

where ^ = ^"^7°, we need to let the following term: 

- 7 • A^)* = J'^A;;,, (4.14) 

be invariant, where 

J^="*7^*, (4.15) 

and 

A'={—^,^). (4.16) 

Since ^'4' is a proper Lorentz scalar and is a proper Lorentz four-current, A'^ must be 
a proper Lorentz four-vector as well. If we define the following correspondence between 
the interaction and external potentials: 

$1 = (*t^)$^ (4.17a) 

A' = (■**)A, (4.176) 

then we see A'^ becomes the Minkowski four-potential A^ = ($, A). 

Considering a four-current eJ^ as a source of four-potential A^, we have a Poisson 
equation 

dpdPA^ = eJ^, (4.18) 
which is the Maxwell equation under a Lorentz gauge 

d^'A^ = 0, (4.19) 

with a current continuity equation derivable from (4.1) 

d^'J^ = 0. (4.20) 

This whole set of the Maxwell equations is invariant under the proper Lorentz transfor- 
mations as well as unitary space-time inversions, although the Minkowski four-potential 
and four-current are "both" time-anticovariant. 
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Take a look at an example: if an electron feels the external force driven by a proton, 
from the correspondence (4.17) and the Maxwell equation (4.18) we get the following 
potential expressions 

^\x) = ^l{x)^e{x) J d^x'[ep¥j,{x')^p{x')]G{x,x'), (4.21a) 

A\x) = '^e{x)^e{x) J d''x'[ep%{x')'y^p{x')]G{x, x'), (4.216) 

where Green's function G{x,x') satisfies 

dpdPG{x,x') = (5W(x - x'). (4.22) 

Hence $^ and A', which depend on not only the wave function of driving source '^p{x') 
but also that of testing body \I'e(a;), are the interaction potentials between the driving 
source and testing body, gaining different meaning from the background potentials $ 
and A driven only by external sources. Though the Minkowski four-potential has been 
successfully used in the Maxwell equation (4.18), the situation in the wave equation (4.1) 
is different: the interaction potentials <I>^ and A' between the spin one- half particles and 
the external electromagnetic fields, should be taken into account, even if <i>^ and A^ do 
not yet form a covariant four-vector. While a time-covariant interaction four-potential 
will be derived in §8. Particularly in the linear approximation, the equation (4.21a) gives 
a point-like static Coulomb potential ee$' = —a/r as we expect. 

From the relations (4.14) through (4.17), the equation (4.1) can also be alternatively 
written down as 

7''P^*(x) = (m + eJ^yl^)*(x). (4.23) 

Here cI'^Af^ happens to be the well-known electro-dynamical interaction potential. So 
unlike the Dirac EM-equation, the constructed equation is a nonlinear wave equation, 
which, in most cases, does not have exact explicit analytical solutions. This result is 
actually a consequence of the Maxwell equation, in which four-potential transforms in the 
same way as four-current under the Lorentz transformations, leading to the appearance 
of a Lorentz invariant electro-dynamical interaction potential in (4.23). 

Originally the Dirac EM-equation was set up to deal with the interaction between 
a point-like electron and an external electromagnetic field, while the internal size of 
electron was ignored. In contrast, the nonlinear equation (4.23) may become valid for 
extended electron, when the interaction potential is expressed by a convolution between 
four-current and four-potential 

V{x) =eJ^*A^ = e j d'^x' J^'{x')Ap{x - x'), (4.24) 

with an integral over the finite 4-d size of extended electron. Using siich a nonsingular 
convolution can avoid the infinity troubles like the infinite self-energy of electron and 
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divergent integrals in Feynman diagrams, caused by the interaction of singular point 
particles. 

If a fermion also participates in other types of interactions, more Lorentz invariant 
interaction terms should be added into V{x). We may just as well generalize (4.23) by 
the following form: 

ij^d^'^ix) = [m + V{x)]'^{x), (4.25) 

where V{x) representing the sum of interaction terms must be invariant under all the 
Lorentz transformations. This equation is basically a nonlinear equation since the in- 
teraction V{x) depends on fermion fields. To completely solve the problem, one needs 
to establish more equations for intermediate bosons by gauge theory. Of course, this is 
only true in principle. When more complicated interactions are involved, solving a com- 
plete set of nonlinear equations are very impractical and much unnecessary, and further 
realizable techniques are required. 

5. ENERGY ORDER OF HYDROGEN-LIKE IONS 

Though (4.1) is a nonlinear equation without exact analytical solutions, we can still 
make a fair approximation that the interaction potentials and are slowly-changing 
functions of particle wavefunctions on account of the small scale of particle itself compared 
with the long-range electromagnetic interaction. Using this approximation and noticing 
the differences between \ipper larger components and lower smaller components of wave 
function, we obtain the non-relativistic Hamiltonian of four major terms after taking 
A' = X r/2 which is valid in the atomic range where B' is approximately constant: 

H.n = l^ + e^^ -i^^-i^x^')-^ir. V$^)(s • 1), (5.1) 

here k = l-l-tr, while 1 = r x p is the orbital angular momentum, s = cr/2 is the spin, and 
the total angular momentum is j = l-|-s. In (5.1), the first term is the kinetic energy; the 
second is the electric potential; the third shows the magnetic Zeeman effect; the forth 
represents the spin-orbit coupling. As we can see, the first three terms are formally the 
same as the conventional results although and A' refer to the interaction potentials 
not the pure external potentials. The forth spin-orbit coupling term has a different sign 
from convention, showing a reversed order of splitting energy levels. 

With the non-relativistic Hamiltonian (5.1), we can write a wave equation idfip = 
HnriP, which is kind of "nonlinear Schrodinger equation" since the interaction poten- 
tials <1>^ and A^ depend on the wave functions. The nonlinear Schrodinger equation has 
certain soliton-like sohitions, and has been applied in many aspects of nonlinear physics. 
However in the derivation of the above non-relativistic Hamiltonian, the positive-energy 
assumption has been taken. This implies: a wave equation with a non-relativistic Hamil- 
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tonian is not invariant under unitary time inversion. In other words, the results ob- 
tained by applying time inversion in a truncated non-relativistic wave equation like the 
Schrodinger equation may not be correct. 

As an example, we can solve the nonlinear relativistic wave equation (4.1) for the 
single-electron model of Hydrogen-like ions under the linear approximation that there is 
a point-like Coulomb potential = —Za/r, but no vector magnetic potential eA^ = 0. 
In this model, (4.1) reduces to 

(7"^"^ + -^)*(a;) = m^'(x), (5.2) 

which differs from the Dirac equation by a matrix /? on the Coulomb term. A more 
general form than (5.2) for scalar central potential (C/r) has been investigated in detail 
(Greiner [15] [16]). Following the standard procedures of solving this type of equations, 
we arrive at the bound state energy solutions of (5.2): 

En, = ±m|l (5.3) 

I [n - (i + 1/2) + V(i + 1/2)2 + (^^)2]2 ; 

which gives a symmetric positive-negative energy spectrum and does not have any diffi- 
culty when Za > 1. 

If Za <^1, we can expand the positive Enj in powers of Za 

{Zaf {ZaY 1 1 
2n2 2rfi + An' 

Where the first term is the rest mass of electron, the second term is the classical quantum 
mechanics result, and the third relativistic term is different from the conventional one 
by a crucial sign change. Although the fine structure splitting energy happens to be the 
same as usual 

A ^ / • • N m(Zar ,1 1 , . ^, 

the energy order for quantum number j with fixed n is changed. Our conclusion is: the 
smaller the total angular momentum j, the larger the positive energy Enj in a certain 
level n of Hydrogen-like ions. For example, the energy level 2Pi/2 is higher than the 
energy level 2P3/2 in the doublet 2P1/2 — 2P3/2 of Hydrogen-like ions, contrary to the 

conventional order. The common practice in experiment is: first detect a spectrum 
of energy levels, then arrange these energy levels in the order derived from the Dirac 
equation. Up till now, the energy order of the doublet of Hydrogen-like ions, has not 
been precisely verified in experiment as far as we know. Let us examine if it is possible 
to double check it by experiment. 

The fine structures of the fines H„(A6563,n = 3^2) and H^(A4861,n = 4^2) 
in the Balmcr series of Hydrogen have been carefiilly investigatc^d (Lewis and Spcdding 
[17]; Spedding et al. [18]), of which no more than two strong components can be resolved 



En, = -[1 - ^ + '-^{-m^ - j::)] + ouzaf]. (5.4) 
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within experimental accuracy. This type of experiments, however, can only show the 
energy difference between two levels, i.e., 2Pi/2 ^ 2P3/2 of Hydrogen atom. It can not, 
by itself, indicate the energy order; which line corresponds to which transition. The fine 
structure of the line A4686(n = 4 ^ 3) of singly ionized Helium He+ has also been studied 
in detail (Paschen [19]). However, if the new energy order is assigned to those energy 
levels, the fit to the experimental data is not worse in consideration of experimental 
accuracies. On the other hand, the Paschen-Back [20] effect by applying strong magnetic 
fields shows a symmetric triplet optical spectrum, and does not make any difference if 
the order of energy levels is reversed. 

On the other hand, the Zeeman effects of multi-electron atoms and ions have been 
vastly investigated (White [21]; Kuhn [22]), which show both "normal order", the bigger 
j the higher energy level, and "abnormal order" , the smaller j the higher energy level. 
When a number of electrons arc involved in multi-electron atoms or ions, more com- 
plicated effects must be taken into account, for example, electrostatic screening, orbit 
penetration and Fermi-Dirac statistics properties of multi-electron systems, and more 
complete equations like the self-consistent Hartree-Fock equations need to be set up. 
They are beyond the scope that the single-electron equation can describe. 

In principle, only the Zeeman effects for Hydrogen-like ions under weak magnetic 
fields can directly and clearly give us the answer of the energy order of Hydrogen-like ions, 
since the asymmetric spectral lines for different transitions may be presented (White [21]; 
Kuhn [22]). The principal fine structure doublet splitting (2Pi/2 — 2P3/2) of Hydrogen 
atom is less than half wave number, hence it is very difficult to detect more splitting 
levels under weak magnetic fields in such a narrow band with good resolution. For this 
reason, the Zeeman effect of Hydrogen atom has not been observed (White [21]), and 
the energy order of the doublet 2P1/2 ~ 2P3/2 in Hydrogen is still an unsolved puzzle 
in experiment. If high resolution (say, one tenth of wave number) in analyzing spectral 
lines is realized, the detection of the Hydrogen Zeeman effect will be easier. For the 
other Hydrogen-like ions such as He+, Li++, Be+++, etc., having bigger fine structure 
splitting energies (oc Z^) as seen in (5.5), the Zeeman effects under weak magnetic fields 
seem detectable, even if most of the simple transitions (n = 2 — > 1, or 3 — > 2) drop into 
ultraviolet or X-ray range, not seen in the range of visible light spectrum. By modern 
techniques in analyzing laser optical spectrum, these experiments are believed realizable. 

From (5.3), the positive ground state energy for n = 1 and j = 1/2 is 

TO 

which approaches zero when Z ^ 00. While the conventional ground state energy of a 
Hydrogen-like ion has a singularity at Za ~ 1: 



Econ = msJX - (Za)^, (5.7) 
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which has no real meaningful solution when Za > 1. The percentage difference between 
(5.6) and (5.7) is 

l-%^ = l-V'l-(^")4, (5.8) 

which is 1% when Z ~ 50, 5% when Z ^ 75 and 15% when Z ~ 100. One can see that 
the errors are significant only for high Z Hydrogcn-likc ions which arc unstable on the 
other hand. This causes real difficulty in doing experiments. An alternative way is to 
measure the ground states of high Z atoms instead of ions, by assuming the outer shell 
electrons have little screening effects on the inner ground state electrons, which may be 
approximately described by single-electron model. So a possible experiment is to use 
laser beam to pump out the electrons of high Z atoms. One may find out the largest 
ionization energies and compare them with the theoretical results (5.6) and (5.7). 

In summary, the energy solutions (5.3) differ from conventional results in a subtle way. 
The differences may not be noticeable if no special attention is paid to it in experiment. 
In fact, the fine structure splitting is exactly the same as usual. Only the energy order 
in a certain level n of Hydrogen-like ions is reversed. This is why we emphasize the 
importance of a thorough experimental investigation of the Zeeman effects of Hydrogen- 
like ions. Ironically no such experiment can be found by far. 

6. PARTICLE CONJUGATION 

After a simple check, we find that charge conjugation defined by = i'j'^^*, does not 
turn (4.23) into the other one with opposite charge. So we need to try something else. 
Now making a complex conjugation * on both sides of (4.23), we get 

-^''*P^^*{x) = (m + eJ^'A^)^*{x). (6.1) 

Considering 7^* = —7^ and 7** = 7*(i = 0, 1, 3), we may pick a unitary transformation 

Os = (6.2) 

which is a complex conjugation operator in spinor space: 

0,Y*0:^ = 7". (6.3) 

Making transformation — Og on both sides of (6.1), we get 

7"^^.(0.*(x)) = (-TO - eJ''^^)(0^*(x)), (6.4) 

where Oo; is a complex conjugation in common space, defined by 

0(^*(a;) = ij°^^-f^^*{x). (6.5) 
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Comparing (6.4) with (4.23), we find that 0^, changes a particle with charge e and mass 
TO into another one with opposite charge — e and opposite mass —to. So we hke to call 
Ou "particle conjugation". 

In §3, we have mentioned that common momentum has two eigenvalues ±to. So we 
may just as well write another equation: 

7'^P^4'_„(x) = -TO*_„(x), (6.6) 

which gives the same type of solutions with the same energy spectrum as the Dirac 
equation (3.5). Adding an interaction term, we get an equation similar to (4.23) but 
with "opposite" mass: 

rP^^_m{x) = {-m + e J'^A^)*_„(x), (6.7) 
which can be transformed under "particle conjugation" into 

j'^P^iO^^.mix)) = (to, - eJ'M^)(0^*_™(a;)). (6.8) 

Clearly (6.8) and (4.23) show opposite charge but same mass. 

Negative charge has been found in experiment long ago, but no negative mass. In 
present-day cosmology, the existence of negative mass is highly controversial. Somehow it 
involves the precise definition of mass. It would be naive to draw any definitive conclusion 
just from this single-electron model. Hence we consider this "particle conjugation" as 
a speculative thought, though we do feel equations (4.23), (6.4), (6.7) and (6.8) are all 
legitimate forms in certain sense. 

Due to our different definitions of unitary time inversion and "particle conjugation", 
the conventional CP or CPT theorem does not hold in our theory. It is possible to find 
more conventional theorems that may not be derivable in our nonlinear theory. We may 
either derive similar ones as replacements, or find completely new theorems. 



7. PERTURBATION TREATMENT 

As long as the interaction potential V{x) is weak enough in (4.25), we may follow Feyn- 
man's perturbation treatment (Feynman [23]), to derive the lowest-order differential cross 
section: 

where subscripts i and / represent the incoming and outgoing waves respectively, u and 
u are the spinors, q = Pf — Pi is the four- momentum transfer and 

V{q) = J d^a;y(x)e*«-^, (7.2) 

is the Fourier amplitude of the interaction potential V{x). 
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In the case of Coulomb interaction between the incident electron beam and nuclein 
target of nuclear number Z: V{x) = — Zq;/|x|, energy is conserved: = Ei. Wo have 
y((j') = —Za/c^ where = (pf — pi)^ = — 4p^(sin|)^. For any unpolarizod incident 
electron beam, the cross section is a sum over final spin states and an average over initial 
spin states. So we get 

da Z^aW 

where the spin sum 

Si = ^ E HPs,s))u{pust)? = + (7-4) 

is different from the conventional result (Bjorken and Drell [7]) since our Coulomb term 
as in (5.2) does not have a matrix /3. Ignoring a constant factor 1/16, finally we have 



2 



(7.6) 



(7.7) 



In the non-relativistic case |p| <C m, it reduces to the Rutherford formula: 

da _ Z'^o?rr? 
dn " p4 (sin 1)4" 

And in the relativistic limit |p| m, it turns out to be 

da _ Z^g^ 
dfl ~ p2 (sin 1)2' 

which differs from the Mott [24] formula by a factor [cot(6'/2)]2. 

Usually the Coulomb interaction is applicable only in low energy scattering, not in 
high energy scattering due to the recoil of nuclear target. In the electron-proton elastic 
scattering experiments (McAllister and Hofstadter [25]), the incident electron energy 188 
Mcv is in the same order as the proton mass 938 Mev. Considering the recoil effect 
of hydrogen, we would better utilize a current-potential interaction with a nonsingular 
convolution (4.24): 

V{x)=eJ^{x)*Al{x). (7.8) 

Its Fourier amplitude becomes 

V{q) =ej d^x'J^ixY^-""' J d^xA^ix - a;')e*«-(^-^') = eJ^{q)Al{q). (7.9) 

From the Maxwell equation (4.18), we may express the four-potential driven by proton 
as follows: 

A^ix) = G{x) * epJP^ix) = J d%G{x - y)epJf:{y), (7.10) 
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where G{x) is Green's function defined by (4.22). In analogy to (7.9), we get 

AP{q) = e,G{q)Jfi{q), (7.11) 
where G{q) = —{q^ + ie)~^. Finally we write down 

V{q) = -e'j^{q)G{q)J^^{q) 



w? M"^ _ -1 _ 

" ~"Y 'eJeIv eJW^^^^' Sf)'-f>'u{pi, Si) ^ i^ <Pf^ Sf)j^u{Pi, Si). (7.12) 

In the relativistic elastic scattering as shown in Fig. 1, the four-momentum transfer 
squared is q^ = ip' — p)"^ = —4EfEi[sm |)^, and the finite mass recoil factor is Ef/Ei = 
[1 + {2E/M)[sm^)^]~^ where E = Ei. For any unpolarized electron beam, the electron 
spin is random over space-time. So when wc calculate the interaction between electron 
and proton, their spins are chosen randomly. The incoming electron spin Si in calculating 
the Feynman diagrams is not related to the incident electron spin s° in the past. Same is 
true for the outgoing electrons. To calculate the total differential cross section per solid 
angle, we sum over all spin states independently: 

da ME^ 

where Si, an extra term that the conventional theory does not have, is given by (7.4): 
El = {EfEi/m^) {sin |)^ when E ^ m, and S2 is the same as in the conventional theory 
(Bjorken and DreU [7]): 

^2 = ^E E \u{Pf,Sf)ru{p^,s,)u{Pf,Sf)^,u{P„Si)\' 



s/,Sj Sf,Si 



^^''-(cos?)^[l-^(tan^)^]. (7.14) 



Combining all these results we obtain 



= ^(cot^)^[ . ], (7.15) 

which differs from the conventional result by a factor [sin(0/2)]~^. 

In comparison with the experimental data (McAllister and Hofstadter [25]) as shown 
in Fig. 2, it is hard to tell which one is better. More advanced research need to be made. 

Generally two form factors need to be considered due to the finite size and anomalous 
magnetic moment of proton (Cahn and Goldhabcr [26]). The modification due to the 
internal structure of proton is significant in high energy scattering processes: the higher 
the incident electron energy, the larger the modification. In the extreme relativistic 
limit ]p] M, the result (7.15) is not valid, and in the certain range jpj M, it is 
not accurate, but in the low relativistic limit m <C |p| ^ M, it may be good enough. 
The lowest-order scattering is quite preliminary. So we would not draw any definitive 
conclusions as yet. 
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8. GAUGE INVARIANCE 



Prom the correspondence (4.17) and the following relation 

"$($1 _ ^ . a')* = J'^yl^* = ■*7''C**)A^*, (8.1) 

we can also write (4.1) or (4.23) as follows: 

7''(P^ - eAl)^ = m*, (8.2) 

by introducing an interaction four-potential 

Al = {^^)A,. (8.3) 

Here A'^ is "time-covariant" , since A^^ is "time-anticovariant" and changes sign 
under unitary time inversion. The equation (8.2) implies a formal similarity to the Dirac 
equation (3.5) if we define an "effective" time-covariant four-momentum 

P^Jf=P,-eAl (8.4) 

or an "effective" time-covariant four-derivative 

d;" = d^+ieAl. (8.5) 

In classical electrodynamics, observable physical quantities can be expressed in terms 
of electromagnetic field strengths: E = — V$ — 9tA and B = V x A which seem invariant 
under gauge transformation for an arbitrary scalar function f{x) (Bjorken and Drell [27]; 
Aithison and Hey [28]): 

A^^A^ + d^f. (8.6) 

However there are several restrictions that need to be taken into account seriously, if 
we want to preserve both gauge invariancc and Lorentz invariance of the Maxwell field 
equation and relativistic fcrmion field equation. 

First of all, as we have pointed out, Minkowski four-potential A^ is time-anticovariant 
(3.12) but four- derivative is apparently time-covariant under time inversion: 

{dt,d^)^{-dt,d^). (8.7) 

So gauge function f(x) must change sign under time inversion to keep the transformed 
four-potential + 9^/ time-anticovariant as same as Af,. Second, gauge function f{x) 
needs to be a proper Lorentz invariant to ensure the transformed four-potential + 
covariant under the proper Lorentz transformations. Third, by making a local phase 
transformation on fermion field 

* *exp[-ia(x)], (8.8) 
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and a local gauge transformation (8.6) on Minkowski field simultaneously in (8.2), we 
arrive at another constraint 

d^a - e(¥*)9^/ = 0. (8.9) 
The general solutions satisfying the above three conditions are 

f{x) = ^ C„(**)", (n = ±1, ±3, ±5, . . .), (8.10a) 

a(a;) =e(V ^^C„ (**)"+! -C7_i log I'^^'l), (n = 1, ±3, ±5, . . .), (8.106) 
^ — ' n + 1 

where C„ are arbitrary real constants to cover the whole phase space. Note: ^ 0, 
otherwise, there is no interaction term in the equation (8.2). 

Furthermore, the Lorentz gauge (4.19) leads to another constraint on the choice of 
gauge function f{x): 

d^d^fix) = 0. (8.11) 

Using (4.23) or more generally (4.25) with (4.24), we can prove step by step that this is 
true if the gauge function takes the form of (8.10a). First we have an expansion 

9^a''(**) = ■*(a^9''*) + (9^9''*)* + 2a^'*a'**. (8.12) 

From (4.25) it comes out 

^i^d^^ = -i{m + V)^. (8.13) 

Making a hermitian conjugation on both sides and then multiplying it by 7'^ from right, 
we get by the relations ^'^^^'^ = ^o^;* ^jj^j ^ _ ^t^o 

^^'^-f^' = i{m + Vy^. (8.14) 

From (8.13) it is derived 

S^S''* = i^a^i'^S^* = -ij'^id^V)^ - (m + Vf'if. (8.15) 
Repeating the procedures in deriving (8.14) leads to 

d^d^"^ = i^j^di^V - (m + y)2*. (8.16) 

Furthermore we have 

25^*5''4' = 2d^'^9^''du^ = 5^*(7''7'' + 7''7'')a^*. (8.17) 
Combining (8.13) and (8.14), we get 

9^*7''7''9^* = (m + Vf^^. (8.18) 

Then we need to find the second term in (8.17). Now let us make a transformation on 
all fermion fields (note 7'^7^t^'^ = 7^) 

=7''*, ='*7'', (8.19) 
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then (8.13) becomes 

Yduil^^) = -i{m + V')!"*- (8-20) 

Here V' is the interaction after the transformation, which for example can be expressed 
by its Fourier amplitude (7.13) in the electron-proton scattering 

T./. N / m2 / M2 f d*q -1 _i„.^ 

xu{pf,Sf)j''Yruipi, Si)u{Pf,Sf)Y^pYu{Pi, Si). (8.21) 

Given = ±1, it is clear to see (8.21) does not change when 7'* and ^ switch twice 
in both electron and proton (or any fermion) currents. So the interaction will not change 
after the transformation (8.19): 

V'{x) = V{x). (8.22) 
Combining (8.20), (8.22) and (8.14), we have 

d^^Yl'^d^^ = -i{m + F)a^*7''* = (m + F)^**. (8.23) 

With (8.18) and (8.23), then (8.17) becomes 

29^*9'** = 2(m + F)^**. (8.24) 

Finally inserting (8.15), (8.16) and (8.24) into (8.12) we arrive at 

a^a'' (■**) = 0, (8.25) 

leading to the constraint (8.11) if the gauge function is given by (8.10a). From the above 
proof, we can see that the gauge condition is quite stringent. It appears that the Lorentz 
gauge is inherent in our theory. Note: the Coulomb gauge, as a special case of the Lorentz 
gauge for static fields, is included. 

In summary, with the explicit solutions (8.10) for f{x) and a{x), the local gauge 
transformations in our nonlinear QED look like: 

* ^ * exp [-ia{x)], (8.26a) 

A^^A^ + d^f{x), (8.266) 

or Al^Al + ^di,a{x), (8.26c) 

which ensure both gauge invariance and Lorentz invariance of the nonlinear equation (8.2) 
and the Maxwell equation with the Lorentz gauge. Note: our nonlinear QED involving 
higher-order fermion fields, is different from the one involving higher-order electromag- 
netic fields, and is also different from the one involving self-interaction potential in the 
Dirac EM-equation, mentioned in some literatures (see reference [29]). 
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We may also construct a nonlinear QED Lagrangian density 

L{x) = ^7^(ia^ - ^eAl)<f - m^* - ^F^^Fi^^m, (8.27) 

here = — d^j^Ai, is the second-rank antisymmetric electromagnetic field ten- 

sor. The Hamiltonian density changes sign under unitary time inversion, so does the 
Lagrangian density. All Euler-Lagrange field equations can be obtained by Hamilton's 
principle. Taking an infinitesimal arbitrary variation on the electromagnetic field 
we get 

d"[Fi,,{li-^)] = e(*7M*)C**), (8.28) 

then keeping the fermion field unchanged d'^{^^) — > 0, we eliminate F^^d'^ {^'^) from 
(8.28) and cancel on both sides to obtain the Maxwell equation d^F^^ = eJ^, that 
reduces to (4.18) under the Lorentz gauge (4.19). Similarly by taking an infinitesimal 
change on the fermion field J^*, we get 

n^^M* - ^67"^;^* - Je(*7''*)^^* - - If^^F^"'^ = 0, (8.29) 

then keeping the electromagnetic field unchanged d^A^^ 0, we eliminate the last term 
in (8.29) to obtain the nonlinear equation (8.2) by the relation (8.1). 

Though the gauge transformation (8.26) does not keep the Lagrangian density (8.27) 

invariant, it docs keep invariant the integral of the Lagrangian density since the extra 
term in the integral makes no contribution: J d^xJ^^d^a = J d'^xd^{J^a) = 0. With the 
definition (8.3), the Maxwell equation (4.18) and the Lorentz gauge (4.19) or (8.25), we 
get 

dpdPA^ = (**) J^^* + 2dp(^'^)dPA^. (8.30) 

Suppose in the free field case (J^^* = 0), the fluctuations of fermion field and photon 
field are small, namely 

dp{^^)=ep{x){^^), (8.31a) 
dPA^ = Kp{x)Ai„ (8.316) 
where €p{x) and kP{x) are small fiuctuation functions. Then we end up with 

dpO^Al + M^(a;)Aj = 0, (8.32a) 

fj,'^{x) = -2ep{x)KP{x). (8.326) 

If i-i^{x) > 0, it is a Klein-Gordon equation for vector boson with a "fluctuation mass" 
H{x), i.e, a harmonic oscillator with a fluctuation energy-momentum. If ii^{x) < 0, it 
shows the decay and oscillation of the interaction field. Both cases are complicated: the 
time-covariant intermediate vector boson field A^, namely the coupled interaction field 
between the Dirac field \1/ and the Minkowski field A^, fiuctuates in space-time. 
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By (8.27) and (8.1), the interaction Hamiltonian density can be written as 




(8.33) 



If the fluctuation is small, we may simply take zeroth-ordcr approximation by setting 



linearize our nonlinear QED to fit experimental data just as accurately as conventional 
linear QED. Furthermore, we may investigate the nonlinear aspects of the theory, by 
considering as a massive intermediate boson, or considering J^A^^ as a nonsingular 
convolution shown in (4.24) to avoid singularities, and rewrite (8.27) as: 



With the basic field quantization techniques, we can then deal with various problems in 
nonlinear QED starting from this Lagrangian density. 

Orthodox QED, on the basis of the Dirac EM-equation, is just a linear theory. In 
this linear QED, electron is point-like and the self-energy of electron blows up when its 
size reduces to zero. Also at each vortex in Feynman diagrams, there is a singular term 
that causes ultraviolet divergences and requires renormalization. In contrast, a nonlinear 
equation of the type (4.25) with (4.24) gives certain "soliton-likc" solutions, not "point- 
like" solutions. So the self-energy of electron would not go to infinity in our nonlinear 
QED. It appears at first sight that this theory is not renormalizable for there is a fourth- 
power term of fermion field in the Lagrangian. However the interaction potential eJ^^^A^ 
in (8.34) is considered as a nonsingular convolution, each vortex in Feynman diagrams 
is smeared, not singular any more as pictured in Fig. 1. Consequently, only two powers 
are left in the perturbation expansion with the other two being integrated out in each 
vortex. Hence this theory is free of ultraviolet divergences and is also renormalizable. 
These are just some general observations. Due to the scope limit of this paper, I leave 
these fundamental issues open for future research. 

So much has been discussed, it is however not my intention to "solve" this nonlinear 
QED in such single paper. Rather my emphasis is on its "derivation" by the implemen- 
tation of unitary time inversion. A complete understanding of nonlinear QED entails 
much more studies both theoretically and experimentally. It has been recognized for a 
long time that nonlinearity may result in fundamentally new phenomena (see reference 
[29] for a history overview). There is a typical problem posed for nonlinear theories: 
since the linear superposition principle of quantum mechanics is not valid in nonlinear 
field equations, it becomes quite problematic to make linear perturbation expansions. A 
great many efforts have been made in the past several decades, to seek non-perturbative 
techniques in solving nonlinear quantum field theories. Many interesting and important 
issues are still far from being settled. 



/x(x) — > 0. In this limit, we may treat Aj^ 



as a time-covariant massless photon and 




(8.34) 
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9. REMARKS 



The theory estabhshed in this paper is strictly Hmited to the Minkowski flat space-time, 
where the Lorentz group is the fundamental symmetry group. Naturally we would ponder 

on the possibility of extending our theory into curved space-time. In a globally curved 
space-time, the Lorentz group is only locally preserved. A conventional treatment on 
the Dirac fields in such a space-time is to use a covariant derivative with spin affine 
connection in the Dirac equation, as outlined in my dissertation (Jin [30]). As far as the 
electromagnetic interaction is concerned, a logical step is to include an electro-dynamical 
interaction potential rather than external potential in the Dirac equation, similar to what 
I have shown here in this paper. After all, a curved space-time with Lorentzian signature 
reduces locally to Minkowski space-time as a limit. Specifically in a static space-time 
with a time-independent metric, the time component of spin afHne connection vanishes, 
we can separate time from space and anticipate a positive-negative symmetric energy 
spectrum (Jin [31] [32]). However in a general space-time with a time-dependent metric, 
time and space are not separable, time inversion is not applicable, and it is no longer 
possible to obtain a positive-negative symmetric energy spectrum. 
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FIG. 1 : Revised Feynman diagram of e-p scattering 
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FIG. 2: Differential cross sections of e-p scattering 



